First-Principles Calculation of Mg(0001) Thin Films: Quantum Size Effect and 

Adsorption of Atomic Hydrogen 
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We have carried out first-principles calculation of Mg(0001) free-standing thin films to study the 
oscillatory quantum size effect exhibited in the surface energy, work function, interlayer relaxation, 
and adsorption energy of the atomic hydrogen adsorbate. The quantum well states have been 
shown. The calculated energetics and interlayer relaxation of clean and H-adsorbed Mg films are 
clearly featured by quantum oscillations as a function of the thickness of the film, with oscillation 
period of about 8 monolayers, consistent with recent experiments. The calculated quantum size 
effect in H adsorption can be verified by observing the dependence of H coverage on the thickness of 
Mg(OOOl) thin films gown on Si(lll) or W(110) substrate which has been experimentally accessible. 

PACS numbers: 73.61.-r, 73.20.At, 73.21.Ac, 
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I. INTRODUCTION 

When the thickness of thin metal films approaches the 
nanoscale, the oscillatory quantum size effects (QSE) as- 
sociated with electronic confinement and interference will 
occur— iSi2il due to the splitting of the energy-level spec- 
trum into subbands normal to the plane of the films. 
Confinement of electrons often leads to strongly modi- 
fied physical properties. It has been shown that a change 
in film thickness by just one atomic layer can result in 
property variations on the order of 1/JV, where N is the 
thickness of the film in terms of monolayers (ML). The 
oscillatory QSE have long been clearly observed in ul- 
trathin metal overlayers on metal substrates^, mostly in- 
volving noble metals. On the other hand, thin metal 
films on semiconductors may be the basis for novel de- 
vices utilizing quantum- well states (QWS). Thus recent 
systematic experimental and theoretical investigation of 
the QSE has mainly been focused to Pb films deposited 
on Si(lll)^a^. u i 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 . 20 i 21 substrate. 

In this paper we present a detailed first-principles 
study of the electronic structure and adsorption ener- 
getics of Mg(OOOl) free standing films. The QSE in some 
other free-standing metal films such as A l 22 ' 23 ' 24 , Li 2 ^, 
and P b 26 i 27 have been theoretically reported in previ- 
ous references. The present study is directly motivated 
by the recent experimental demonstration that highly 
perfect ultrathin epitaxial Mg(OOOl) films can be grown 
on Si(lll) substrate by low-temperature deposition and 
annealing 2 ^, and on W(110) substrat o 29 ' 30 ^. Us- 
ing angle-resolved photoemission spectroscopy (ARPES) 
technique, Aballe et ali^ have showed that each time 
a QWS state falls in the wave-vector and energy range 
of the upper branch of the Mg sp band, a new peak is 
visible in the photoemission spectrum, with the thick- 
ness interval between the two sequential peaks being 
about 8 ML. Schiller et al— have extensively measured 
the electronic structure of magnesium from Mg(OOOl) 
monolayers to bulk. More recently, Koitzsch et al. and 
Schiller et al. have reported the spin-splitting effects 
in ultrathin Mg(OOOl) films due to the coupling of the 



Mg(0001)/W(110) interface electronic structure and the 
QWS states. It is expected that further work related 
with the Mg(OOOl) thin films will be reported afterwards. 
From this aspect a thorough theoretical investigation of 
the QSE in Mg(OOOl), including the energetics and the 
interlayer relaxation, is necessary and will be helpful for 
the future experimental reference. It should be men- 
tioned that the theoretical study of QSE on Mg(OOOl) 
was initiated by Feibelman^ 2 .. 

The other object in this paper is to study the QSE 
character in the atomic adsorption energetics, which has 
been neglected by most of the previous studies. Since 
the adsorption property is closely characterized by the 
chemical bonding between the adsorbate and the surface 
of the substrate, thus when the substrate is ultra-thin, 
the QSE in the substrate will also influence the behavior 
of the surface adsorption. Here as a case study we choose 
the atomic hydrogen as the adsorbate on Mg(OOOl), since 
H is the most simple element, also since the influence 
of atomic hydrogen adsorption on the surface electronic 
structure of the metals has been extensively studied with- 
out emphasis on QSE. Our results show that in the ul- 
trathin Mg(OOOl) films, the adsorption energy of atomic 
H displays a well-defined QSE. 

This paper is organized as follows: In Sec. II, the ab 
initio based method and computational details is out- 
lined. In Sec. Ill, the surface properties of the Mg(OOOl) 
films, including the electronic structure, surface energy, 
work function, and interlayer relaxation, as a function of 
the thickness of the films, are presented and discussed. 
Also the properties of adsorption of atomic hydrogen 
monolayer onto Mg(OOOl) surface is discussed in detail 
by presenting the sensitivity of the adsorption energy to 
the thickness of the Mg(OOOl) films. Finally, Sec. IV 
contains a summary of the work and our conclusion. 



II. COMPUTATIONAL METHOD 

The calculations were carried out using the Vi- 
enna ab initio simulation packag based on density- 
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FIG. 1: GGA energy bands and density of electron states 
(right panel) of hep bulk Mg. The dashed line denotes Fermi 
level. 

functional theory with ultrasoft pseudopotentials 3 - 3 - and 
plane waves. In the present film calculations, free- 
standing Mg(0001) films in periodic slab geometries were 
employed. The periodic slabs are separated by a vac- 
uum region equal to 20 A. In all the calculations be- 
low, a surface (lxl) was employed for the supercell 
slab. The Brillouin-zone integration was performed us- 
ing Monkhorse-Pack scheme 34 with allxllxlfc- 
point grid, and the plane-wave energy cutoff was set 250 
eV. Furthermore, the generalized gradient approximation 
(GGA) with PW-91 exchange-correlation potential has 
been employed with all atomic configurations fully re- 
laxed. First the total energy of the bulk hep Mg was 
calculated to obtain the bulk lattice constants. The cal- 
culated a- and c-lattice parameters are 3.201 Aand 5.186 
A, comparable with experimental values of 3.21 Aand 
5.213 A, respectively. The use of larger fc-point meshes 
did not alter these values significantly. A Fermi broaden- 
ing of 0.1 eV was chosen to smear the occupation of the 
bands around Ep by the Fermi-Dirac function. 

III. RESULTS AND DISCUSSION 

Band structure. — We first studied the properties of 
electronic structures of Mg(0001) films. As a first step, 
we present in Fig. 1 the band structure and the density of 
states (DOS) of bulk hep Mg. One can see that the DOS 
of bulk Mg is nearly free-electron like (~ y/e). This is 
different from its close neighbor Be whose DOS resembles 
somewhat that of a semiconductor due to the absence of 
core electrons in Be atoms. Also one can see From Fig. 1 
that there are two filled state at T with energies around 
1 eV, while in the case of Be the corresponding states are 
above the Fermi energy. Although the outmost electronic 
configuration of elemental Mg is 3s 2 , one can see from 
Fig. 1 that the p-orbital component in bulk Mg plays as 
well an important role around Ep. 

The electronic structure properties of Mg(0001) film 
are shown in Fig. 2 as a function of the thickness of 
the film. Here we only plot the energy levels at Y point 
without interlayer relaxation. The energy zero is set at 




Thickness N (ML) 

FIG. 2: Calculated (GGA) energies at V in Mg(0001) thin 
films as a function of thickness, with the energy set to zero 
at the Fermi level. The right panel replot the bulk energy 
dispersion in the [0001] direction. 

the Fermi level of each film. Interlayer relaxation effect 
has also been studied and it is found that the overall 
thickness dependence of the energies is similar to that 
without relaxation. For comparison and illustration, the 
bulk energy dispersion along the [0001] (T — » A) direction 
[see Fig. 1], which determines the energy range for the 
quantum well states, is plotted again in Fig. 2 (right 
panel). One can see that the energy gap at T is small 
due to the hybridization of Mg s and p orbitals, with 
the energy of ~0.5 eV. The QWS states arise from the 
upper band. When the thickness of the film is increased 
to ~8 ML, then a QWS state, with the energy crossing 
the Fermi level, occurs. The next energy crossing with 
the Fermi level occurs at the film thickness of 16 ML. Our 
calculated results of QWS states are in good agreement 
with recent experimental ARPES measurement. 

For sp metals the QWS states are often analyzed in the 
framework of the phase accumulation mode l 36 ! 37 . Here 
the free-standing Mg(0001) film is considered as a quan- 
tum well confining electrons between the two vacuums 
in the slab. Only such kj_ (perpendicular component of 
bulk wave vector) values of the electrons are allowed that 
fulfill the stationary state condition for integer n, 

2k ± Nd + 2$ = 2im, (1) 

where N is the number of atomic layers in the film, 
d = c/2 the interlayer spacing, and $ the phase shift 
of the electronic wave function upon reflection at the 
film- vacuum interface. Using Eq.(l) one can calculate 
the periodicity for the QWS states crossing the Fermi 
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FIG. 4: Work function of clean and H-adsorbed Mg(0001) 
thin films as a function of thickness. 



FIG. 3: (a) Monolayer energy E(N)/N, (b) corresponding 
energy difference AE(N), and (c) surface energy for fully re- 
laxed Mg(0001) 1 x 1 slabs as a function of thickness. 



level, AN = 7r/(fc{d), where is the perpendicular 
component of Fermi wave vector. From the right panel 
in Fig. 2, one can see that the upper branch of the bulk 
sp band runs through about 25% of the Brillouin zone, 
k{ = 0.25tt/c. One gets AiV = 8. Therefore a new 
QWS state occurs every 8 ML, which is verified in the 
left panel in Fig. 2 that an energy branch moves down, 
crossing the Fermi level for every incremental increase in 
the film thickness of 8 ML. 

Energetics. — Figure 3(a) shows the total energy 
E t (N)/N per ML as a function of the thickness of the 
Mg(0001) film. The atoms in the slabs have been fully 
relaxed during calculations. One can see from Fig. 3(a) 
that with increasing the thickness, E t (N)/N gradually 
approaches a constant value which in the limit is equal 
to the energy per atom in the bulk Mg. 

An energetic quantity more suitably tailored to the 
QSE is the surface energy which is defined as one-half of 
the energy difference between the film and the bulk with 
the same number of atoms, including the proper subtrac- 
tion of a term linear in TV 38 . The thickness dependence 
of surface energy is shown in Fig. 3(b). It reveals that 
consistent with the result of the electronic structure in 
Fig. 2, the surface energy follows a simple oscillatory 
form with the period of A7V ~ 8. These oscillations arise 
from the occupation of electronic levels close to Fermi 
surface. Also one can see that the oscillation pattern 
is not as good as one expects. This is due to the fact 
that the atomic arrangement in the slabs has been fully 
relaxed during the calculation. 

Figure 4 (line with triangles) shows the work function 
as a function of the thickness of Mg(0001) film for relaxed 
atomic geometry. One can see the work function is also 
featured by an oscillatory form multiplied by a damping 
factor. The oscillation period in thickness is again AN ~ 
8. 



TABLE I: Interlayer relaxations given in percent, Ad^i+i, of 
Mg(0001) as a function of the thickness of the film. 
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Interlayer relaxation. — The interlayer relaxation, 
Adj^+i, is given in percent with respect to the unrelaxed 
interlayer spacing, do, i-e., Adj^+i = 100(di.i+i — d^)/d,Q. 
d^i+i is the interlayer distance between two adjacent lay- 
ers parallel to the surface calculated by total energy min- 
imization, do = c/2 is the bulk interlayer distance. As 
mentioned above, all layers in the slab were allowed to 
relax. Obviously, the signs + and — of Adi^+i indicate 
expansion and contraction of the interlayer spacings, re- 
spectively. The relaxation of Mg(0001) surface as a func- 
tion of the thickness of the film is summarized in Table I. 
Furthermore, the interlayer relaxations are also plotted 
in Fig. 5 for more clear illustration. One can see: (i) 
The two outmost layers relax significantly from the bulk 
value, in agreement with experimental observation^. In 
the whole range of layers that we considered, the topmost 
interlayer relaxation is always outward (Adi^ > 0). The 
value of Adi^ starts from +8% for a slab with only 2 
ML, and approaches a final value of with increasing 
the thickness of Mg(0001) film. In contrast to the behav- 
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FIG. 5: Interlayer relaxations of Mg(0001) thin films as a 
function of thickness. 




FIG. 6: The four different adsorption sites for H adatom on 
Mg(000I) surface. 



ior of AcZi t 2, the second interlayer relaxation is always 
inward (Ad2,3 < 0). Note that the first interlayer sep- 
aration on most metal surfaces is contracted, Mg(0001) 
is one of the few exceptions; (ii) The interlayer spacings 
oscillate as a function of the thickness of the film with 
a damped magnitude. The oscillation period is about 8 
ML in the thickness, thus clearly indicating the QSE in 
the interlayer relaxation. After 30 ML, which is the max- 
imal layers considered here, the oscillations are invisible, 
which suggests that the semi-infinite surface limit is now 
reached. 

Adsorption of atomic hydrogen: QSE of binding 
energy. — To further illustrate the physical properties in- 
fluenced by finite size of the thin films, in this section we 
focus our attention to the adsorption of atomic hydrogen 
on Mg(0001) thin films. To the best of our knowledge, 
the reflection of QSE by the adsorption features has been 
neglected by most of the previous studies. 

Before we study the hydrogen adsorption properties as 
a function of the thickness of the film, we need to deter- 
mine the energetically favorable adsorption site. Since 
the preference of adsorption site is not sensitive to the 
thickness of the substrate, thus to look for this prefer- 
ence, it is sufficient to give a study on the slabs with 
fixed thickness of Mg(0001) substrate, which at present 
is chosen to be 9 ML. We choose four most probable 
adsorption sites, namely, on-top (Tl), bridge (B2), fee 
(H3), and hep (T4) sites, which are schematically in- 
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FIG. 7: Calculated (a) binding energy of H adatom and (b) 
adsorbate height as a function thickness of Mg(0001) films. 



dicated in Fig. 6. The binding energy is calculated 
using the following equation: Binding energy [atomic 
H]= -(£[H/Mg(0001)] - £[Mg(0001)] - 2£[H])/2 where 
£'[H/Mg(0001)] is the total energy of a slab which con- 
sists of 9 layers of Mg atoms and one H atom on each side 
keeping inversion symmetry, £7[Mg(0001)] is total energy 
of the slab without H atoms, and E\H] is total energy of a 
free H atom which is put in a supercell with the size of 10 
A. As a result, the calculated binding energy of atomic H 
for different adsorption configurations is 1.317 eV (Tl), 
2.394 eV (H3), and 2.319 eV (T4). The B2 site is unstable 
and after the zero-temperature relaxation, the H adatom 
at B2 will diffuse to H3 site. Thus the fee (H3) site is 
most stable and in the following discussions, the atomic 
H is always put on H3 site during the simulation. Note 
that this H3 preference for H adsorption on Mg(0001) is 
different from the cases in Be(0001) and W(100) wherein 
the bridge-site is preferre d 41 i 42 . 

After finding the preferred atomic H adsorption site 
(H3), we gave a series of calculations for the binding en- 
ergy of the H adsorbate as a function of the thickness 
of Mg(0001) thin films. The results are summarized in 
Fig. 7(a). One can see that the binding energy curves 
up at small film thickness, followed by damped oscilla- 
tions when increasing the Mg ML in the slab. Thus the 
binding energy of atomic H depends on the thickness of 
the quantum films in an oscillatory way. The oscilla- 
tion period in thickness is about 8 monolayers, indicat- 
ing a well-defined QSE in the adsorption of atomic H on 
Mg(0001). In experiment this QSE of atomic adsorption 
can be observed by investigating the dependence of H 
coverage on the monolayers of Mg(0001) thin films. Also 
we have calculated the H adsorbate height and the results 
are plotted in Fig. 7(b), which again shows the periodic 
oscillations indicative of QSE. Furthermore, we have cal- 
culated the work function for the H-adsorbed Mg(0001), 
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which is shown in Fig. 4 (line with squares). One can see 
that compared to the clean Mg(0001), the work function 
is reduced by the presence of H adlayer, implying that the 
charge is transferred from H to Mg. Also it shows in Fig. 
4 that as in the case of clean surface, the work function 
of H-adsorbed Mg(0001) is oscillatory in the amplitude 
with respect to the thickness of the film. The oscillation 
period in thickness is again AN ~ 8. 

IV. CONCLUSION 

In summary, the Mg(0001) thin films have been studied 
by density-functional theory pseudopotential plane-wave 
calculations. The dependence of electronic structure, en- 
ergetics, and interlayer relaxation upon the thickness of 



the film has been fully investigated, clearly showing the 
metallic QSE of the film, in consistent with the recent 
experiments. We have also studied the atomic hydrogen 
adsorption on the Mg(0001) film. It has been shown that 
the adsorption energy of H oscillates with the increase of 
Mg monolayers in the slab. As the other energetic quan- 
tities, this oscillation in the adsorption energy can also 
be explained by the occurrence of the QWS states in the 
ultrathin metal films. 



Acknowledgments 

This work was supported by the CNSF No. 10544004 
and 10604010. 



1 J.J. Paggel, T. Miller, and T.-C. Chiang, Science 283, 1709 
(1999). 

2 Z. Tesanovic, M.V. Jaric, and S. Maekawa, Phys. Rev. 
Lett. 57, 2760 (1986). 

3 A. E. Meyerovich and S. Stepaniants, Phys. Rev. Lett. 73 
316 (1994). 

4 A. Kawabata, J. Phys. Soc. Jpn. 62, 3988 (1993). 

5 T. Miller, A. Samsavar, G.E. Franklin, and T.C. Chiang, 
Phys. Rev. Lett. 61, 1404 (1988); H. Zeng and G. Vidali, 
ibid. 74, 582 (1995); B. J. Hinch, C. Koziol, J. P. Toennies, 
and G. Zhang, Europhys. Lett. 10, 341 (1989). 

6 I. B. Altfeder, D.M. Chen, and KA. Matveev, Phys. Rev. 
Lett. 80, 4895 (1998). 

7 V. Yeh, L. Berbil-Bautista, C.Z. Wang, K.M. Ho, and M.C. 
Tringides, Phys. Rev. Lett. 85, 5158 (2000). 

8 W.B. Su, S.H. Chang, W.B. Jian, C.S. Chang, L.J. Chen, 
and T.T. Tsong, Phys. Rev. Lett. 86, 5116 (2001). 

9 LB. Altfeder, V. Narayanamurti, and D.M. Chen, Phys. 
Rev. Lett. 88, 206801 (2002). 

10 P. Czoschke, H. Hong, L. Basile, and T.-C. Chiang, Phys. 
Rev. Lett. 93, 036103 (2004). 

11 M.H. Upton, CM. Wei, M.Y. Chou, T. Miller, and T.-C. 
Chiang, Phys. Rev. Lett. 93, 026802 (2004). 

12 DA. Ricci, T. Miller, and T.-C. Chiang, Phys. Rev. Lett. 
93, 136801 (2004). 

13 J.H. Dil, J.W. Kim, S. Gokhale, M. Tallarida, and K. Horn, 
Phys. Rev. B 70, 045405 (2004). 

14 Y. Guo et al, Science 306, 1915 (2004). 

15 M.H. Upton, T. Miller, and T.-C. Chiang, Phys. Rev. B 
71, 033403 (2005). 

16 P. Czoschke, H. Hong, L. Basile, and T.-C. Chiang, Phys. 
Rev. B 72, 035305 (2005). 

17 P. Czoschke, H. Hong, L. Basile, and T.-C. Chiang, Phys. 
Rev. B 72, 075402 (2005). 

18 Y.-F. Zhang, J.-F. Jia, T.-Z. Han, Z. Tang, Q.-T. Shen, 
Y. Guo, Z. Q. Qiu, and Q.-K. Xue, Phys. Rev. Lett. 95, 
096802 (2005). 

19 Z.Y. Zhang, Q. Niu, and C. K. Shih, Phys. Rev. Lett. 80, 
5381 (1998). 

20 Juarez LF. Da Silva, Phys. Rev. B 71, 195416 (2005). 

21 E. Ogando, N. Zabala, E.V. Chulkov, and M.J. Puska, 
Phys. Rev. B 71, 205401 (2005). 



22 P.J. Feibelman, Phys. Rev. B 27, 1991 (1983). 

23 A. Kiejna, J. Peisert, and P. Schraroch, Surf. Sci. 432, 54 
(1999). 

24 S. Ciraci and LP. Batra, Phys. Rev. B 33, 4294 (1985); LP. 
Batra, S. Ciraci, G.P. Srivastava, J.S. Nelson, and C.Y. 
Fong, Phys. Rev. B 34, 8246 (1986). 

25 J.C. Boettger and S.B. Trickey, Phys. Rev. B 45, 1363 
(1992); U. Birkenheuer, J.C. Boettger, and N. Rosch, 
Chem. Phys. Lett. 341, 103 (1995); K.F. Wojciechowski 
and H. Bogdanow, Surf. Sci. 397, 53 (1998); K. Doll, N.M. 
Harrison, and V.R. Saunders, J. Phys.: Condens. Matter 
11, 5007 (1999). 

26 CM. Wei and M.Y. Chou, Phys. Rev. B 66, 233408 (2002). 

27 P. Saalfrank, Surf. Sci. 274, 449 (1992). 

28 L. Aballe, C. Rogero, and K. Horn, Phys. Rev. B 65, 
125319 (2002). 

29 F. Schiller, M. Heber, V.D. Servedio, and C. Laubschat, 
Phys. Rev. B 70, 125106 (2004). 

30 C. Koitzsch, C. Battaglia, F. Clerc, L. Despont, M.G. Gar- 
nier, and P. Aebil, Phys. Rev. Lett. 95, 126401 (2005). 

31 F. Schiller, R. Keyling, E.V. Chulkov, and J.E. Ortega, 
Phys. Rev. Lett. 95, 126402 (2005). 

32 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

33 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

34 H.J. Monkhorst and J.D. Pack, Phys. Rev. B 13, 5188 
(1976). 

35 V.M. Amonenko, V. Ye. Ivanov, G.F. Tikhinskij, and VA. 
Finkel, Phys. Met. Metallogr. 14, 47 (1962). 

36 P.M. Echenique and J.B. Pendry, J. Phys. C 11, 2065 
(1978). 

37 N.V. Smith, Phys. Rev. B 32, 3549 (1985). 

38 J.C. Boettger, J.R. Smith, U. Birkenheuer, N. Rosch, S.B. 
Trickey, J.R. Sabin, and S.P. Apell, J. Phys.: Condens. 
Matter 10, 893 (1998), and references therein. 

39 F.K. Schulte, Surf. Sci. 55, 427 (1976). 

40 Ismail and E.W. Plummer, Bull. Am. Phys. Soc. 43, 1, 46 
(1998). 

41 R. Stumpf and P.J. Feibelman, Phys. Rev. B 51, 13748 
(1995). 

42 M.R. Barnes and R.F. Willis, Phys. Rev. Lett. 41, 1729 
(1978). 



